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Abstract. We present a new fully dynamic algorithm for maintaining betweenness centrality (BC) 
of vertices in a directed graph G = {V, E) with positive edge weights. BC is a widely used parameter 
in the analysis of large complex networks. We achieve an amortized ■ log^ n) time per update, 

where n = \V\ and v* bounds the number of distinct edges that lie on shortest paths through any single 
vertex. This result improves on the amortized bound for fully dynamic BC in [21122] by a logarithmic 
factor. Our algorithm uses new data structures and techniques that are extensions of the method in 
the fully dynamic algorithm in Thorup [2S] for APSP in graphs with unique shortest paths. For graphs 
with V* = 0{n), our algorithm matches the fully dynamic APSP bound in Thorup |28| . which holds 
for graphs with v* = n — 1, since it assumes unique shortest paths. 


1 Introduction 

Betweenness centrality (BC) is a widely-used measure in the analysis of large complex networks, and 
is defined as follows. Given a directed graph G = {y,E) with \V\ = n, \E\ = m and positive edge 
weights, let axy denote the number of shortest paths (SPs) from x to y in G, and Oxyiv) the number 
of SPs from x to y in G that pass through v, for each pair x, y € R. Then, BC{v) = t^v 

As in [S], we assume positive edge weights to avoid the case when cycles of 0 weight are present in 
the graph. 

The measure BC{v) is often used as an index that determines the relative importance of v in 
G, and is computed for all v £ V. Some applications of BC include analyzing social interaction 
networks m, identifying lethality in biological networks [20] , and identifying key actors in terrorist 
networks msi . In the static case, the widely used algorithm by Brandes |5| runs in 0{mn+v? logn) 
on weighted graphs. Several approximation algorithms are available: mM for static computation 
and, recently, m for dynamic computation. Heuristics for dynamic betweenness centrality with 
good experimental performance are given in [10116126] . but none provably improve on Brandes. 
The only earlier exact dynamic BC algorithms that provably improve on Brandes on some classes 
of graphs are the recent separate incremental and decrementatj algorithms in HSUS]. Recently, 
we give in |2I] (see also [22|) the first fully dynamic algorithm for BC (the PR method) that is 
provably faster than Brandes for the class of dense graphs (where m is close to n^) with succinct 
single-source SP dags. Table [1] contains a summary of these results. 

In this paper, we present an improved algorithm for computing fully dynamic exact betweenness 
centrality: our algorithm FFD improves over PR by a logarithmic factor using data structures and 
technique that are considerably more complicated. Our method is a generalization of Thorup [28] 
(the Thorup method) which computes fully dynamic APSP for graphs with a unique SP for every 
vertex pair; however, a key step in BC algorithms is computig all SPs for each pair of vertices {all 
pairs aM shortest paths — APASP). We develop a faster fully dynamic algorithm for APASP, which 
in turn leads to a faster fully dynamic BC algorithm than PR. 

* This work was supported in part by NSF grant CCF-1320675. Authors’ affiliation: Dept, of Computer Science, 
University of Texas, Austin, TX 78712. 

^ Incremental/decremental refer to the insertion/deletion of a vertex or edge; the corresponding weight changes that 
apply are weight decreases/increases, respectively. 














Our Results. Let i'* be the maximum number of distinct edges that lie on shortest paths through 
any given vertex in G] for convenience we assume i/* = I7(n). 

Theorem 1. Let U be a sequence of L2{n) fully dynamic vertex updates on a directed n-node graph 
G = (F, E) with positive edge weights. Let u* bound the number of distinct edges that lie on shortest 
paths through any single vertex in any of the updated graphs or their vertex induced subgraphs. 
Then, all BC scores (and APASP) can be maintained in amortized time ■ log^ n) per update 

with algorithm ffd. 

Similar to the PR algorithm, our new algorithm FFD is provably faster than Brandes on dense 
graphs with succinct single-source SP dags. It also matches the Thorup bound for APSP when 
u* = 0{n). 

Our techniques rely on recomputing BC scores using certain data structures related to shortest 
paths extensions (see Section [2|) . These are generalizations of structures introduced by Demetrescu 
and Italiano [8] for fully dynamic APSP (the DI method) and Thorup, where only one SP is 
maintained for each pair of vertices. Our generalizations build on the tuple-system introduced in [19] 
(the NPRdec method) for decremental APASP (see Section [3|), which is a method to succinctly 
represent all of the multiple SPs for every pair of vertices. Our algorithm also builds on the fully 
dynamic APASP (and BC) algorithm in PR, which runs in amortized dog^ n) cost per update 

(see Section [3]). Finally, one of the main challenges we address in our current result is to generalize 
the ‘level graphs’ of Thorup to the case when different SPs for a given vertex pair can be distributed 
across multiple levels. 


Paper 

Year 

Time 

Weights 

Update Type 

DR/UN 

Result 

Brandes 

2001 

0(mn) 

NO 

Static Alg. 

Both 

Exact 

Brandes 

2001 

0{mn + 'n? log n) 

YES 

Static Alg. 

Both 

Exact 

Geisberger et al. 0 

2007 

Heuristic 

YES 

Static Alg. 

Both 

Approx. 

Riondato et al. |24| 

2014 

depends on e 

YES 

Static Alg. 

Both 

e-Approx. 


Semi Dynamic 


Green et al. |10| 

2012 

0(mn) 

NO 

Edge Inc. 

Both 

Exact 

Kas et al. |12| 

2013 

Heuristic 

YES 

Edge Inc. 

Both 

Exact 

NPR [n] 

2014 

0(u* -n) 

YES 

Vertex Inc. 

Both 

Exact 

NPRdec 

2014 

0{v*^ ■ logn) 

YES 

Vertex Dec. 

Both 

Exact 

Bergamini et al. [4] 

2015 

depends on e 

YES 

Batch (edges) Inc. 

Both 

e-Approx. 


Fully Dynamic 


Lee et al. 1161 

Singh et al. |26j 
KourtellisH- |14| 
Bergamini et al. 

PR I21I22I 

This paper (ffd) 

2012 

2013 

2014 

2015 
2015 

2015 

Heuristic 

Heuristic 
0(mn) 
depends on e 
■ log3 n) 
0(v*'^ ■ log^n) 

NO 

NO 

NO 

YES 

YES 

YES 

Edge Update 

Vertex Update 
Edge Update 
Batch (edges) 
Vertex Update 
Vertex Update 

UN 

UN 

Both 

UN 

Both 

Both 

Exact 

Exact 
Exact 
e-Approx. 
Exact 

Exact 

Table 1. Related results (DR stands for Directed and Ub 

for Undirected) 


Discussion of the parameters m* and iz*. Let m* be the number of distinct edges in G that 
lie on shortest paths; n*, defined above, is the maximum number of distinct edges on any single 
source SP dag. Clearly, u* < m* < m. 

- m* vs m: In manv cases, m* <C m: as noted in m, in a complete graph (m = 0(n^)) where 
edge weights are chosen from a large class of probability distributions, m* = 0(n log n) with high 
probability. 

- iz* vs m* : Clearly, u* = 0{n) in any graph with only a constant number of SPs between every 
pair of vertices. These graphs are called fe-geodetic [23| (when at most k SPs exists between two 
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nodes), and are well studied in graph theory |27l2ll7j . In fact u* = 0{n) even in some graphs that 
have an exponential number of SPs between some pairs of vertices. 

In contrast, m* can be 0(n^) even in some graphs with unique SPs, for example the complete 
unweighted graph Kn- 

Another type of graph with ly* m* is one with large clusters of nodes (e.g., as described 
by the planted i-partition model [325]). Consider a graph H with k clusters of size n/k (for some 
constant k > 1) with 5 < w{e) < 25, for some constant <5 > 0, for each edge e in a cluster; between 
the clusters is a sparse interconnect. Then m* = I7(n^) but v* = 0(n). (The connections between 
clusters can be arbitrary, thus BC scores are non-trivial to compute.) For all of the above classes of 
graphs, our BC algorithm will run in amortized O(n^) time per update (O hides polylog factors). 
More generally we have: 

Theorem 2. Let U be a sequence of L2{n) updates on graphs with 0{n) distinct edges on shortest 
paths through any single vertex in any vertex-induced subgraph. Then, all BC scores (and APASP) 
can he maintained in amortized time 0{n^ ■ log^ n) per update. 

Our algorithm uses 0{m ■ v*) space, extending the 0{mn) result in DI for APSP. Brandes uses 
only linear space, but all known dynamic algorithms require at least l7(n^) space. 

Overview of the Paper. In Section [2] we describe how to obtain a fully dynamic BC algorithm 
using a fully dynamic APASP algorithm. The remaining sections in the paper are devoted to 
developing our improved fully dynamic APASP algorithm. In Section [3] we review the NPRdec, PR 
and Thorup methods. In Section 0] we introduce the level tuple-system framework for APASP, with 
particular reference to the new data structures specifically developed for our result. In Section [5] 
we present our algorithm FFD, and we describe its main components in detail. Section [6] describes 
two important challenges that arise when generalizing Thorup to APASP setting; addressing these 
two challenges is crucial to the correctness and efficiency of our algorithm. Finally, in Section [7] we 
establish the amortized time bound of ■ log^ n) for ffd and its correctness. 

2 Fully Dynamic Betweenness Centrality 

The static Brandes algorithm [5] computes BC scores in a two phase process. The first phase 
computes the SP out-dag for every source through n applications of Dijkstra’s algorithm. The 
second phase uses an ‘accumulation’ technique that computes all BC scores using these SP dags in 
0(n • n*) time. 

In our fully dynamic algorithm, we will leave the second phase unchanged. For the first phase, 
we will use the approach in the incremental BC algorithm in [18], which maintains the SP dags using 
a very simple and efficient incremental algorithm. For decremental and fully dynamic updates the 
method is more involved, and dynamic APASP is at the heart of maintaining the SP dags. Neither 
the decremental nor our new fully dynamic APASP algorithms maintain the SP dags explicitly, 
instead they maintain data structures to update a collection of tuples (see Section m- We now 
describe a very simple method to construct the SP dags from these data structures; this step is 
not addressed in the decremental NPRdec algorithm and it is only sketched in the fully dynamic PR 
algorithm. 

For every vertex pair x, y, the following sets R*{x, y), L*{x, y) (introduced in DI) are maintained 
in NPRdec, and in both of our fully dynamic algorithms; 

- R*{x, y) contains all nodes y' such that every shortest path x y m. G can be extended with the 
edge {y,y') to generate another shortest path x y ^ y'. 

- L*{x, y) contains all nodes x' such that every shortest path x y m G can be extended with the 
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edge {x\x) to generate another shortest path x' ^ x y. 

These sets allow us to construct the SP dag for each source s using the following simple algorithm 
BUILD-DAG. 


Algorithm 1 build-dag(G, s, w, D) (w is the weight function; D is the distance matrix) 
1: for each t £V do 
2: for each u € R*[s,t) do 

3: if D{s, t) + w(f, u) = D{s, u) then add the edge (t, u) to dag(s) 


In our fully dynamic algorithm R* and L* will be supersets of the exact collections of nodes 
defined above, but the check in Step [3] will ensure that only the correct SP dag edges are included. 
The combined sizes of these R* and L* sets is 0(n • u* • logn) in our ffd algorithm, hence the 
amortized time bound for the overall fully dynamic BC algorithm is dominated by the time bound 
for fully dynamic APASP. In the rest of this paper, we will present our fully dynamic APASP 
algorithm FFD. 

3 Background 

In this section we review prior work upon which we build our results. For each, we highlight the 
inherited notation we use and the main ideas we extend. 

3.1 The NPRdec Decremental APASP Algorithm jl9j 

NPRdec generalizes the decremental APSP algorithm in DI [8] to obtain a decremental algorithm 
for APASP and BC. For the decremental APSP algorithm DI develops a novel method to maintain 
(unique) shortest paths [8]. DI defines an LSP as a path where every proper subpath is a shortest 
path in the graph. By efficiently maintaining all LSPs after each update, DI presents an efficient 
decremental APSP, which is then extended to a fully dynamic APSP algorithm with additional 
tools. This provides a fully dynamic algorithm for APSP that runs in O(n^logn) amortized time 
per update. NPRdec extends this result to APASP by introducing the tuple-system to replace the 
need to maintain every SP and LSP in DI. We now briefly review this system, referring the reader 
to [l9] for more details. Let d{x, y) denote the shortest path length from x to y. 

A tuple, T = {xa, by), represents a set of paths in G, all with the same weight, and all of which 
use the same first edge {x,a) and the same last edge {b,y). If the paths in r are LSPs, then r is 
an LST (locally shortest tuple), and the weight of every path in r is w(x, a) + d{a, b) + 'w{b, y). In 
addition, if d{x, y) = w(x, a) -\- d{a, b) + w( 6 , y), then r is a shortest path tuple (ST). 

A triple 7 = {T,wt, count), represents the tuple r = {xa,by) that contains count number of 
paths from x to y, each with weight wt. NPRdec uses triples to succinctly store all LSPs and SPs 
for each vertex pair in G. 

Left Tuple and Right Tuple. A left tuple (or Atuple), = {xa,y), represents the set of LSPs from 
X to y, all of which use the same first edge {x,a). A right tuple (r-tuple) = {x,by) is defined 
analogously. For a shortest path r-tuple = {x, by), L{Tr) is the set of vertices which can be used 
as pre-extensions to create LSTs in G, and for a shortest path £-tuple T£ = {xa,y), R{t£) is the set 
of vertices which can be used as post-extensions to create LSTs in G. 

NPRdec maintains several other sets such as P* and P for each vertex pair. Since our algorithm 
also maintains generalizations of these sets, we will discuss them when we present the data structures 
used by our algorithm in Section HI 
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Similar to DI, the NPRdec algorithm first deletes all the paths containing the updated node 
using a procedure cleanup, and then updates the tuple-system to maintain all shortest paths in 
the graph using procedure fixup. The main difference, in terms of data structures, is the use of 
tuples to collect paths that share the same first edge, last edge and weight. 

3.2 The PR Fully Dynamic APASP Algorithm [21j 

In our recent fully dynamic PR algorithm for APASP [2T] , we build on the tuple-system introduced 
in NPRdec. PR also incorporates several elements in the DI method of extending their decremental 
APSP algorithm to fully dynamic (though some key elements in PR are significantly different from 

DI). 

One difference between NPRdec and PR is the introduction of HTs and LHTs; these are exten¬ 
sion of historical and locally historical paths in DI to tuples and triples defined as follows (THTs 
and TLHTs are not specifically used in this paper, however we include their definitions for com¬ 
pleteness): 

HT, THT, LHT, and TLHT. Let r be a tuple in the tuple-system at time t. Let t' <t denote 
the most recent step at which a vertex on a path in r was updated. Then r is a historical tuple 
(HT) at time t if r was an ST-tuple at least once in the interval [t',t]', r is a true HT (THT) at 
time t if it is not an ST in the current graph. A tuple r is a locally historical tuple (LHT) at time t 
if either it only contains a single vertex or every proper sub-path in it is an HT at time t; a tuple 
T is a true LHT (TLHT) at time t if it is not an LST in the current graph. 

Similar to DI for unique SPs, PR forms LHTs and TLHTs in its fully-fixup procedure (which 
adapts the fixup procedure in NPRdec to the fully dynamic case) by combining compatible pairs 
of HTs. An important method introduced in DI for efficiency in fully dynamic APSP is the notion 
of a ‘dummy update’ sequence. Extending this method to an efficient algorithm for APASP does 
not appear to be feasible, so in PR, a new dummy update sequence (that uses elements in a different 
method by Thorup [28]) is introduced. PR then defines the Prior Deletion Graphs (PDGs) (that 
are related to the level graphs maintained in Thorup - see below) to study the complexity of the PR 
algorithm. In our algorithm FFD we will use graphs similar to the PDGs in PR; these are described 
in Section IH 


3.3 The Thorup Fully Dynamic APSP Algorithm [ 28 j 

In [28], Thorup improves by a logarithmic factor over DI (for unique shortest paths) by using a level 
system of decremental only graphs. The shortest paths and locally shortest paths are generated 
level by level leading to a different complexity analysis from DI. When a node is removed from the 
current graph, it is also removed from every older level graph that contains it. The implementation 
of the Thorup APSP algorithm is not fully specified in [28]. For our ffd algorithm, we present 
generalizations of the data structures sketched in Thorup (see Section 0] for a summary of these 
data structures). 

4 Data Structures for Algorithm ffd 

Our algorithm FFD requires several data structures. Some of these are already present in NPRdec, PR 
and Thorup, while others are newly defined or generalized from earlier ones. We will use components 
from PR such as the abstract representation of the level system using PDGs (see Section [4]) and 
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the flag bit /3 for a triple, the Marked-Tuples scheme introduced in NPRdec (see for more 

details), and the maintenance of level graphs from Thorup. 

In this section, we describe all data structures used by our algorithm. In Table [2] we summarize 
the structures we use, including those inherited from [19121128] . The new components we introduce 
in this paper to achieve efficiency for fully dynamic APASP, are described in section 14.21 and listed 
in Table O Part D. 


4.1 A Level System for Centered Tuples 

Algorithm FFD uses the PDGs defined in PR as real data structures similar to Thorup for APSP. 
This is done in order to generate a smaller superset of LSTs than PR, and this is the key to achieving 
the improved efficiency. Here we describe the level system and the data structures we use in ffd, 
with special attention to the new elements we introduce. 

As in [8121128] we build up the tuple-system for the initial n-node graph G = (P, E) with n 
insert updates (starting with the empty graph), and we then perform n updates according to the 
update sequence E. After 2n updates, we reset all data structures and start afresh. 

Our level system is a generalization of Thorup to fully dynamic APASP. For an update at step 
t, let k be the position of the least significant bit with value 1 in the binary representation of t. 
Then the t-th update activates the level k, and deactivates all levels j < k hy folding these levels 
into level k. These levels are considered implicitly in PR, and using the same notation, we will say 
that time{k) = t, and level{t) = k] moreover Gt indicates the graph after the t-th update. Note 
that the largest level created before we start afresh is r = log 2 n. 

Centering vertices and tuples/triples. As in Thorup, each vertex v is centered in level k = 
level{t), where t is the most recent step in which v was updated. A path p in a tuple is centered 
in level k' = level{t'), where t' is the most recent step in which p entered the tuple system (within 
some tuple) or was modified by a vertex update. Hence, in contrast to Thorup, a triple can represent 
paths centered in different levels. Thus, for a triple 7 = {{xa, by), wt, count) we maintain an array 
G^ where 


G^\i] = number of paths represented by 7 that are all centered in level i 

It follows that = count. The level center of the triple 7 is the smallest (i.e., most recent 

level) i such that G^\i\ ^ 0 . 

Level graphs (PDGs). The PR algorithm defines PDGs as follow: Let t be the current update 
step, let t' < t, and let W be the set of vertices that are updated in the interval of steps [t' -|- l,t]. 
The prior deletion graph (PDG) Et' is the induced subgraph of Gt' on the vertex set V{Gt') — W. 

In PR, the PDGs are used only in the analysis, and are not maintained by the algorithm. Here, 
in FFD, we will maintain a set of local data structures for each PDG that is relevant to the current 
graph; also, in a small change of notation, we will denote a level graph for time t' < t as Ei^i, 
where k' = level{t') rather than the PR notation of Ft'. These graphs are similar to the level graphs 
in Thorup. As in Thorup, only certain information for T), is explicitly maintained in its local data 
structures: the STs centered in level k plus all the extensions that can generate STs in T),. The data 
structures used by our algorithm to maintain triples are Global and Local, which we now describe. 

Global Structures. The global data structures are P*, P, L and R (see Table EJ Part A). 
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Notation 

Data Structure Appears 

Part A :: Global Data Structures (for each pair of nodes {x,y)) 

P{x,y) 

all (centered) LHT s from x to y with weight as key 

(|H] for paths, [19J for 
LSTs) 

P*{x,y) 

all (centered) HT s from x to y with weight as key 

L{x,by) 

{x' : (x'x, by) denotes a (centered) LHT} 

[21], ([E] for LSTs) 

R{xa,y) 

{y' : {xa,yy') denotes a (centered) LHT} 

m], (Ug for LSTs) 

Marked-Tuples 

global dictionary for Marking scheme 

[2T|T9] 

Part B :: Local Data Structures (for each active level for each pair of nodes (x,y)) 

P:{x,y) 

STs from x to y centered in level i 

(sketched in [28] 
for paths) 

Li{x,y) 

{x' : {x'x, y) denotes an t-tuple for SPs centered in level i} 

Ri{x,y) 

{y' : {x,yy') denotes an r-tuple for SPs centered in level i} 

LC:{x,y) 

the subset {x' G L* (x, y) : x' is centered in level i} 

RCt{x,y) 

the subset {y' G Rpx,y) : y' is centered in level i\ 

dicti 

dictionary of pointers from local STs to global P and P* 

new 

Part C :: Inherited Data Structures 

Pi-y) 

flag bit for the (centered) triple 7 

m 

level (t) 

level activated during t-th update 

m 

time{k) 

most recent update in which level k is activated 

m 

N 

nodes (centered in levels) deactivated in the current step 

m 

Pk 

level graph (PDG) created during time{k)-t\i update 

( 1211 as PDG, [28J for paths) 

Part D :: New Data Structures 


distribution of all paths in triple 7 among active levels 

new 

DL{x,y) 

linked-list containing the history of distances from x to y 

new 

LN{x, y, wt) 

the set {b : 3 {xa, by) of weight wt in P{x, y)} 

new 

RN{x, y, wt) 

the set {a : 3 {xa, by) of weight wt in P{x,y)} 

new 


Table 2. Notation summary 


— The structures P*{x,y) and P{x,y) will contain HTs (including all STs) and LHTs (including 
all LSTs), respectively, from x to y. They are priority queues with the weights of the triple and 
a flag bit /3 as key. For a triple 7 in P, the flag bit /^(y) = 0 if the triple 7 is in P but not in 
P*, and P{'y) = 1 if the triple 7 is in P and P*. 

— The structure L{x,by) {R{xa,y)) is the set of all left (right) extension vertices that generate a 
centered LHT in the tuple-system. 

Local Structures. The local data structures we introduce in this paper are L*,R*, LC* and RC* 
(see Table [21 Part B). These are generalization of the data structures sketched in Thorup for unique 
SPs in the graph. For every pair of nodes {x,y): 

— The structure P*{x^y) contains the set of STs from x io y centered in Pj. It is implemented as 
a set. 

— The structure L*{x,y) {R*{x,y)) contains all left (right) extensions that generate a shortest 
£-tuple (r-tuple) centered in level i. It is implemented as a balanced search tree. 

— The structure LC*{x,y) {RC*{x,y)) contains left (right) extensions centered in level i that 
generate a shortest £-tuple (r-tuple) centered in level i. It is implemented as a balanced search 
tree. 

— A dictionary dicU, contains STs in P* using the key [x,y,a,b] and two pointers stored along 
with each ST. The two pointers refer to the location in P(x, y) and P*{x, y) of the triple of the 
form {x,a,b,y) contained in P*{x,y). 
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In order to recompute BC scores (see Section[2]) we will consider R*{x, y) = |J- R*{x, y) and similarly 
L*{x,y) = \JiL*{x,y). 

4.2 New Structures 

We introduce two completely new data structures which are essential to achieve efficiency for our 
FFD algorithm. Both are needed to address the Partial Extension Problem (PEP) which does not 
appear in Thorup (see section [ 6 j) . 

Distance History Matrix. The distance history matrix is a matrix DL where each entry is a 
pointer to a linked list: for each x,y € V, the linked list DL{x, y) contains the sequence of different 
pairs {wt,k), where each one represents an SP weight wt from x to y, along with the most recent 
level k in which the weight wt was the shortest distance from x to y in the graph Rk- The pair with 
weight wt in DL{x, y) is double-linked to every triple from x to y with weight wt in the system. 
Precisely, when a new triple 7 from x to y of weight wt is inserted in the algorithm, a link is formed 
between 7 and the pair {wt,k) in DL{x,y). With this structure, the ffd algorithm can quickly 
check if there are still triples of a specific weight in the tuple-system, especially for example when 
we need to remove a given weight from DL (e.g. Step [21]- Alg. [5|). Note that the size of each linked 
list is 0 (log n). 

Historical Extension (HE) Sets RN and LN. Another important type of structure we intro¬ 
duce are the sets RN and LN. These structures are crucial to select efficiently the set of restored 
historical tuples that need to be extended (see Section IQ]) . RN{x,y,wt) {LN{x,y,wt) works sym¬ 
metrically) contains all nodes b such that there exists at least one tuple of the form (xx, 6 y) and 
weight wt in P{x,y). Similarly to DL, every time a new triple 7 of this form is inserted in the 
tuple-system, a double link is created between 7 and the occurrence of b in RN{x,y,wt) in order 
to quickly access the triple when needed (e.g. Steps [TH to flGl - Alg. [9|). 

The total space used by DL, RN and LN is 0(n^ log n). This is dominated by the overall space 
used by the algorithm to maintain all the triples in the tuple-system across all levels (see Lemma 
m Section [7|). 

5 The FFD Algorithm 

Algorithm FFD is similar to the fully dynamic algorithm in PR and its overall description is given 
in Algorithm |2| The main difference is the introduction of the notion of levels as described in 
Section (U and their activation/deactivation as in Thorup. At the beginning of the t-th update 
(with k = level(t)), we first activate the new level k and we perform ff-update (Alg. |3|) on the 
updated node v. As in PR and shown in Table[2|- Part C, the set W consists of all vertices centered 
at these lower deactivated levels. All vertices in Af are re-centered at level k during the t-th update 
(Alg. 121 Step[5|), and ‘dummy’ update operations are performed on each of these vertices. Note that 
M contains the 2 ^ — 1 most recently updated vertices in reverse order of update time (from the 
most recent to the oldest). Procedure ff-update is invoked with the parameter k representing 
the newly activated level. Finally, all levels j < k are deactivated (Alg. [21 Step [ 6 ]). 

FF-UPDATE. As in [msEi], the update of a node occurs in a sequence of two steps: a cleanup phase 
and a fixup phase. Here, we call this update ff-update and it is a sequence of 2 calls: ff-cleanup 
and FF-FIXUP. Briefly, ff-cleanup removes all LHPs in the tuple-system containing the updated 
node (see Section [5T]) . while ff-fixup identifies and adds the STs and LSTs in the updated graph 




Algorithm 2 ffd(G, I?, w',/c) 

1: activate the new level k 
2: FF-UPDATE(i;, w', fc) 

3: generate the set Af 

4: for each u £ J\f in decreasing order of update time do 
5: ff-update(u, w', /c) {dummy updates} 

6: deactivate all levels lower than k 


that are not yet in the tuple-system (see Section [5.2p . Both ff-cleanup and ff-fixup are more 
involved algorithms than their counterparts in PR and the resulting algorithm will save a O(logn) 
factor over the amortized cost in PR. 


Algorithm 3 ff-update(u, w', k) 

1: FF-CLEANUP(li) 

2: FF-Fixup(i;, w', fc) 


Dummy Updates. Calling ff-update only on the updated nodes gives us a correct fully dynamic 
APASP algorithm. However, the number of LHTs generated could be very large making this strategy 
not efficient in general. Thus, as in PR, the ffd algorithm performs a sequence of dummy updates 
as follows. Consider an update on v at time t. Let k = level{t). As shown in algorithm [2l M is the 
set of 2^ — 1 most recently updated nodes vt-i,vt- 2 ,' " where j = 2^ — 1. These are the nodes 
centered at levels smaller (more recent) than k before the f-th update is applied. For each vertex u in 
A/", starting from the most recent to the oldest, ffd calls ff-update on the node u. This procedure 
removes all LHTs containing u in every active level using ff-cleanup, and immediately reinserts 
them in the newly activated level k by performing ff-fixup and using the local data structures. 
These are the dummy updates. Dummy updates have the effect of removing from the tuple-system 
any path vr containing a vertex in Af that is no longer an LSP in the current graph, because vr is 
removed by the corresponding ff-cleanup step and will not be restored during ff-fixup. 

We will establish in Section [7] that ffd correctly updates the data structures with the amortized 
bound given in Theorem [TJ 

5.1 Description of ff-cleanup 

FF-CLEANUP removes all the LHPs going through the updated vertex v from all the global structures 
P, P*, L and R, and from all local structures in any active level graph Pj that contains these triples. 
This involves decrementing the count of some triples or removing them completely (when all the 
paths in the triple go through v). The algorithm also updates local dictionaries and the DL, RN and 
LN structures. Algorithm ff-cleanup is a natural extension of the NPRdec cleanup. An extension 
of the NPRdec cleanup is used also in PR but in a different way. 


Algorithm 4 FF-CLEANUP(u) 

1: He ■<— 0; Marked-Tuples 0 

2 : 7 ■<— [(u, v), 0 , 1 ]; C-y [center(v)] = 1 ; add [ 7 , C-y] to He 
3: while He ^ 0 do 

4: extract in S all the triples with the same min-key [wt, x, y] from He 

5: FF-CLEANUP-£-extend(5,[u;t, a:, y]) (see Algorithm [Sj 

6 : FF-CLEANUP-r-extend{5',[iHt, x, y]) 
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FF-CLEANUP starts as in the NPRdec algorithm. We add the updated node v to He (Step [2] 
- Alg. 0]) and we start extracting all the triples with same min-key (Step 0]- Alg. 0]). The main 
differences from NPRdec start after we call Algorithm [5j As in [19], we start by forming a new triple 
7 ' to be deleted (Steps [5]- Alg. 0|). A new feature in Algorithm [5] is to accumulate the paths that we 
need to remove level by level using the array C. This is inspired by Thorup where unique SPs are 
maintained in each level. However, our algorithm maintains multiples paths spread across different 
levels using the arrays associated to LSTs, and the technique used to update the arrays is 
significantly different and more involved than the one described in Thorup. Step [ 6 ] - Alg. [5] calls 
FF-CLEANUP-CENTERS (Alg. [6]) that will perform this task. 


Algorithm 5 FF-CLEANUP-Aextend(S', [wt,x,y\) 

1: for every b such that (xx^by) G S do 

2: let Sf) Q S be the set of all triples of the form (xx, by) 

3: let f count' be the sum of all the counts of all triples in Sb 

4: for every x' in L(x,by) s.t. {x'x,by) ^ Marked-Tuples do 

5: wt' •<— ict -h w{x', x); 7' ■<— {{x'x, by), wt', f count') 

6: Cy ^ FF-CLEANUP-CENTERS(7^ 

7: add Yl',C^i] to He 

8: remove 7' in P{x' ,y) j j decrements count by f count' 

9: set new center for 7" = {(x'x,by),wt') in P(x',y) as argmin^(C^//[2] 7^ 0) 

10: if a triple for (x'x, by) exists in P(x',y) then 

11: insert (x'x, by) in Marked-Tuples 

12: else 

13: delete x' from L(x, by) and delete y from R(x'x, b) 

14: if no triple for ((x' — ,by),wt') exists in P(x',y) then 

15: remove b from RN(x',y,wt') 

16: if no triple for ((x'x,—y),wt') exists in P(x',y) then 

17: remove x trom LN(x',y,wt') 

18: if a triple for ((x'x,by),wt') exists in P*(x',y) then 

19: remove 7' in P*(x',u) jj decrements count by fcount' 

20: if 7^ ^P*(x',y) then 

21: remove the element with weight wt' from DL(x',y) if not linked to other tuples in P*(x',y) 

22: for each i do 

23: decrement [2] paths from 7' G P*(x',y) 

24: if 7' is removed from P*(x',y) then 

25: if x' is centered in level i then 

26: if y j > i, Pj (x,y) = (h then 

27: remove x' from L*(x,y) and remove x' from LC*(x,y) 

28: else if P*(x,y) = 0 then 

29: remove x' from L*(x,y) 

30: if y is centered in level i then 

31: if Vj > 2, Pj(x', fo) = 0 then 

32: remove y from R* (x', b) and remove y from RC* (x', b) 

33: else if P*(x',b) = 0 then 

34: remove y from R*(x',b) 


FF-CLEANUP-CENTERS takes as input the generated triple 'j' of the form [x'x, by) which contains 
all the paths going through the updated node v to be removed, and the set of triples Sb of the form 
{xx,by) that are extended to x' to generate 7 '. This procedure has two tasks: (1) generating the 
Cy vector for the triple 7 ' that will be reinserted in He for further extensions, and ( 2 ) updating the 
Cyi vector for the tuple 7 " in P{x',y) (note that 7 " is the corresponding triple in P of 7 ', before 
we subtract all the paths represented by 7 ' level by level). 

(1) - This task, which is more complex than the second task (which is a single step in the algorithm, 
see point (2) below), is accomplished in steps 0] to El Alg. [ 6 ] and uses the following technique. In 
step 0] - Alg. [ 6 l we store into the log n-size array C the distribution over the active levels for the 
set of triples in Sb that generates 7 ' using the left extension to x'. In order to generate the correct 
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Algorithm 6 FF-CLEANUP-CENTERS(7', Sb) 

1 : let y = ({x'x, by), wt', /count') (the triple of the form {{x'x, by), wt') that contains all the paths through v to be removed) 
2: let 7 " = ({x'x, by), wt', /count") (the triple of the form {{x'x, by), wt') in P{x', y). Note that 7 ' represents a subset of 7 ") 
S: j argmaxj {C^f/ [j] ^ 0) j j This is the oldest level in which a path in 7 " appeared for the first time 

4: C' ^ ~ 1; • • ■ ) 0] // This is the sum (level by level) of the triples of the form {xai, by) that go through v and 

are extending to 7 ' during this stage 
5: create a new center vector C^/ for the triple 7 ' as follows 
6 : for all the levels m > j we set C^/ [m] = 0 

7: for the level j we set C^/ [j] = ^ ^'[^] 

8 : for all the levels i < j we set C^f [i] = C'[i] 

9: [r — 1, . .., 0] C^ff [r — 1, ..., 0] — [r — 1, .. ., 0] // We update the C vector for 'y" G P{x', y) 

10: return C^/ // We return the correct vector for the generated 7 ' triples 


vector Cy (to associate with the triple 7 '), we need to reshape the distribution in C according to 
the corresponding distribution of the triple 7 " € P. The reshaping procedure works as follows; we 
first identify the oldest level j in which the triple 7 " appeared in P for the first time (Step [3]- Alg. 
ED- Recall that we want to remove 7 ' paths containing v from 7 ", and 7 " does not exist in any level 
older than j. Vector C is the sum of for all 'y G Sb (StepE]- Alg. [ 6 |). Those triples are of the 
form {xai,by) and they could exist in levels older, equal or more recent than j. But the triples in 
Sb that were present in a level older than j, were extended to 'y" in P for the first time in level j. 
For this reason, step [7]- Alg. [ 6 ] aggregates all the counts in C in levels older or equal j in C^f[j]. 
Moreover, for each level i < j, if a triple 7 G S';, is present in the level graph Tj with count paths 
centered in level i, then Pi also contains its extension to x' that is a subtriple of 'y" located in level 
i with at least count paths. Thus for each level i < j step [ 8 ] - Alg. El copies the number of paths 
level-wise. This procedure allows us to precisely remove the LHPs only from the level graphs where 
they exist. After C is reshaped into Cy (steps El to El - Alg. El), the algorithm returns this correct 
array for 7 ' to Alg. El 

( 2 ) - This task is performed by the simple step El Alg. El which is a subtraction level by level of 
LHPs. 

After adding the new triple j' to He (Step [7]- Alg. E]), the algorithm continues as the NPRdec 
(Steps [7] to EHl - Alg. E]) with some differences: we need to update centers, local data structures, 
DL, RN and LN. We update the center of 7 ' using Cy (Step El- Alg. E|)- If is a shortest triple, 
we decrement the count of 7 ' G P*{x',y) (Step EHl - Alg. E])- If 'y' is completely removed from 
P*{x',y) and DL(x',y,wt') is not linked to any other tuple, we remove the entry with weight wt' 
from DL{x',y) (Step EH- Alg. E])- Moreover, we subtract the correct number of paths from each 
level using the (previously built) array Cy (Step [23]- Alg. E])- Finally for each active level i, if 7 ' is 
removed from P*{x',y), we take care of the sets L* and R* i Steps (241 to (Ml - Alg. E])- In the process, 
we also update LC* and RC* in case the endpoints of 7 ' are centered in level i. If 7 ' is completely 
removed from P{x', y), using the double links to the node b in RN(x', y, wt'), we check if there are 
other triples that use b in P{x',y) (Step EH- Alg. El): if not we remove b from RN{x',y,wt'). A 
similar step handles LN{x',y,wt'). 


5.2 Description of ff-fixup 

FF-FIXUP is an extension of the fixup in PR rather than NPRdec. This is because of the presence 
of the control bit (3 (defined in Section E]) , and the need to process historical triples (that are not 
present in NPRdec). Algorithm FF-FIXUP will efficiently maintain exactly the LSTs and STs for each 
level graph in the tuple-system. This is in contrast to PR, which can maintain LHTs that are not 
LSTs in any level graph (PDG). ff-fixup maintains a heap Hf of candidate LHTs to be processed 
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in min-weight order. The main phase (Alg. [7|) is very similar to the fixup in PR. The differences are 
again related to levels, centers and the new data structures. 

We start describing Algorithm!?! We initialize Hf hy inserting the edges incident on the updated 
vertex v with their updated weights (Steps [ 2 | to [?|- Alg. [ 8 ]), as well as a candidate min-weight triple 
from P for each pair of nodes (x, y) (Step[T0]- Alg.[ 8 |). Then we process Hj hy repeatedly extracting 
collections of triples of the same min-weight for a given pair of nodes, until Hf is empty (Steps [3| 
todOj- Alg. [7|). We will establish that the first set of triples for each pair {x,y) always represents 
the shortest path distance from x to y (see Lemma ED, and the triple extracted are added to the 
tuple-system if not already there (see Alg. [9| and Lemma [ 8 D . As in PRfor efficiency, among all the 
triples present in the tuple-system for a pair of nodes, we select only the ones that need to be 
extended: this task is performed by Algorithm [9] (this step is explained later in the description). 
After the triples in S are left and right extended by Algorithm 1101 we set the bit = 1 for 

each triple 7 ' that is identified as shortest in S, since 7 ' is correctly updated both in P*{x,y) and 
P{x,y) (Stepllj- Alg. [7|). Finally, we update the DL{x,y) structure by inserting (or updating if an 
element with weight wt is already present) the element with weight wt and the current level at the 
end of the list (Step [TO] - Alg. [7|). This concludes the description of Algorithm [7| 

We now describe Algorithm [9| which is responsible to select only the triples that have valid 
extensions that will generate LHTs in the current graph. In Algorithm [9l we distinguish two cases. 
When the set of extracted triples from x to y contains at least one path not containing v (StepEj- 
Alg. ED, then we process all the triples from P{x,y) of the same weight. Otherwise, if all the paths 
extracted go through v (Step [HI- Alg. ED, we only use the triples extracted from Hj. 


Algorithm 7 FF-Fixup(u, w', /c) 

1: Hf ■<— 0; Marked-Tuples 0 

2: FF-POPULATE-HEAP(u, w', fc) 

3: while Hf^% do 

4: extract in S' all the triples with min-key [wt^x^y] from Hf 

5: if 5' is the first extracted set from Hf for x,y then 

6: S ff-new-paths(5', P(x, y)) 

7: FF-FixuP-^-extend(5,[u;t, X, y]) (see Algorithm 

8: FF-FixuP-r-extend(S',[ir’t,a:,y]) 

9: for every 7 € 5 set ^('y) = 1 

10: add an element with weight wt and level k to DL(x,y) or update the level in the existing one 


Algorithm 8 FF-POPULATE-HEAp(u, w', k) 

1 : for each (u,u) do 
2: w(u, u) = w'(u, u) 

3: if w(u,u) < oo then 

4: 7 = {{uv, uu), w{'U, v), 1 ); [k] ■<— 1 

5: update-num( 7 ) ■<— curr-update-num; num-v-paths( 7 ) •<— 1 

6 : add [ 7 , 07 ] to Hf and P{u,v) 

7: add u to L(—, uu) and u to —) 

8 : for each (u, u) do 

9: symmetric processing as Steps [2}{3 above 

10 : for each x,y £V do 

11: add a min-key triple [ 7 , C7] £ P{x,y) to Hf 


Both cases have a similar approach but here we focus on the former which is more involved than 
the latter. As soon as we identify a new triple 7 ' we compute its center j by using its associated 
array Cy (Step [4]- Alg. ED. This is straightforward if compared to ff-cleanup where we first need 
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Algorithm 9 FF-NEW-PATHS(5', Pxy) 

1: 5 •<— 0; let i be the min-weight level associated with DL(x,y) 

2: if P*{x,y) increased min-weight after cleanup then 
3: for each 7 ' € 5 with-key do 

4: let y = {{xa',b'y),'wt, count') and j = argminj(C^/[j] 7 ^ 0) 

5: if 7 ' is not in P*{x,y) then 

6 : add 7 ' in P*{x,y) and S; add x to L*{a',y) and y to R*(x,b') 

7: add b' to RN{x^y,wt)] place a double link between 7 ' and DL{x^y^wt) 

8: else if 7 ' is in P{x^y) and P*(x^y) with different counts then 

9: replace the count of 7 ' in P*{x, y) with count' and add 7 ' to S 

10: add 7 ^ to Pj(x,y) and dictj 

11: add x to Lj(a',y) and y to Rj(x,b') 

12: add x to LCj{a',y) {y to RC^ {x,b')) if x (y) is a level i center 

13: add 7 ^ in S 

14: for each b' E RN{x,y,wt) do 

15: if 3h < i : LJ^{x,b') ^ (ii then 

16: add any 7 ' of the form (xx^b'y) and weight wt in P*{x,y) with ^{'y') = 1 to 5 

17: else 

18: for each 7 ' E S' containing a path through v do 

19: let 7 ' = {{xa', b'y), wt, count') and k the current level 

20: add ■y' with paths( 7 ^,i;) to P*{x,y), and [ 7 ^,C-y/] to 5" 

21: add 7 ' to P^{x,y) and dictk, x to L'^{a',y) and y to R'^{x,b') 

22: add x to LC^{a',y) (y to RC'^{x,b')) if x {y) is a level k center 

23: return 5" 


to update the center arrays. We add this triple to P*{x,y) and to S, which contains the set of 
triples that need to be extended. We also add 7 ' to Pj{x,y) (Steps [ 10 ] and [ 21 ]- Alg. [9]). We update 
dictj to keep track of the locations of the triple in the global structures. A similar sequence of steps 
takes place when all the extracted paths go through v (Steps [T 8 ]to [ 22 ]- Alg. [9|). The only difference 
is that the local data structures to be updated are only the T^ data stuctures (Steps [ 21 ] and [ 22 ] - 
Alg. E]). 

A crucial difference from PR and this algorithm is the way we collect the set S of triples to be 
extended. Here we require the new HE data structures RN and LN (see Section 14. 2|) because of 
PEP instances (see Section [ 6 ]). Let i be the min-weight level associated with DL(x,y). Eor each 
node b E RN{x, y, wt) we check if L'^ix, b) contains at least one extension, for every h < i lSteps[T41 
to[in]- Alg. [9]). In fact we need to discover all tuples with /? = 1 that are inside a PEP instance. In 
this instance, the triples restored as STs may or may not be extended. We cannot afford to look at 
all of them, thus our solution should check only the triples with an available extension. Moreover, 
all the extendable triples with with /3 = 1 have extension only in levels younger than the level 
where they last appear as STs. Thus, we check for extensions only in the levels h < i. 

Using the HE sets, is the key to avoid an otherwise long search of all the valid extensions for 
the set of examined triples with /3 = 1. In particular, without the HE sets, the algorithm could 
waste time by searching for extensions that are not even in the tuple-system. Correctness of this 
method is proven in section [7] After the algorithm collects the set S of triples that can be extended, 
FF-FIXUP calls FF-Fixup-Uextend lAlg. [T0|1 . 

Here we describe the details of algorithm 1101 Its goal is to generate LHTs for the current graph 
G by extending HTs. Let h be the center of Sf, defined as the most recent center among all the 
triples in Sb, and let j be the level associated to the first weight wt' larger than wt in DL{x,y). 
The extension phase for triples is different from PR: in fact, the set of triples Sb could contain only 
triples with /^(q) = 1. In PR, the corresponding set Sb contains only triples with /^(q) = 0. We 
address two cases: 

(a) - If Sb contains at least one triple q with /3(q) = 0, we extend Sb using the sets L* and R* with 
h <i < j (Steps [7] to [16]- Alg fTO]) . In fact, the set Sb contains at least one new path that was not 
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Algorithm 10 FF-FixuP-^-extend(5,[i(;t, x, y]) 

1: for every b such that (xx^by) G S do 

2: let Sb Q S be the set of all triples of the form (xx, by) 

3: let fcount' be the sum of all the counts of all triples in Sb‘, let h be the center(Sb) 

4: if 37 E 5b : ^{'y) = 0 then 

5: let j be the level associated to the minweight wt' > wt in DL(x,y) 

6 : for every active level h < i < j do 

7: for every x' in L*{x,b) do 

8: if {x'x,by) ^ Marked-Tuples then 

9: wt' <— wt + w(x', x); 7 ' <— {(x'x, by), wt', f count') 

10: Cy <— FF-FIXUP-CENTERS(5b); add 7 ' to Hf 

11: if a triple 'y" for {{x'x,by),wt') exists in P{x',y) then 

12: update the count of 7 " in P{x', y) and C^// = -h C^/ 

13: add {x'x,by) to Marked-Tuples 

14: else 

15: add [ 7 ^C'^/] to P{x',y); add x' to L{x,by) and y to R{x'x,b) 

16: set I3{'y') = 0 ; set update-num{ 7 ') 

17: for every level i < h do 

18: for every x' in LC*{x,b) do 

19: execute steps fSl to 1161 

20 : else 

21: let j be the level associated to the minweight wt in DL{x,y) 

22 : for every level i < j do 

23: for every x' in LC*{x, b) do 

24: execute steps to 1161 


Algorithm 11 FF-FIXUP-CENTERS(5;,) 

1: let C' = C'y be the sum (level by level) of the new paths that are found shortest 

2 : let j be argmaxj(C'[j] 7 ^ 0 ), and k = center{x') 

3: if fc < j then 

4: for all the levels i < k we set C^/ [i] = C'[i] 

5: for the level k we set C^/ [fc] = 

6 : for all the levels m > k we set Cy [m] = 0 

7: else 

8 : = C' 

9: return C^/ 
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extended in the previous iterations when wt was the shortest distance from x to y (because of the 
f3{'y) = 0 triple). The LST generated in this way remains centered in level h. Moreover we extend 
Sb also using the sets LC* and RC* with i < h (Steps [T7] to [19]- Alg JTOl) . This ensures that every 
LST generated in a level i lower than h is centered in i thanks to the extension node itself. This 
technique guarantees that each LHT generated by Algorithm [101 is an LST centered in a unique 
level. 

(b) - In the case when there is no triple 7 in Sb with f3{'y) = 0, then there is at least one extension 
to perform for Sb and it must be in some level younger than the level where wt stopped to be the 
shortest distance from x to y (this follows from the use of the HE sets in Alg. [9]). To perform these 
extensions we set j as the level associated with the min-weight element in DL{x, y), and we extend 
Sb using the sets LC* and RC* with i < j (Steps Ell to [Ml- Alg JTO]) . Again, every LHT generated 
is an LST centered in a unique level. Finally, following PR, every generated LHT is added to P and 
Hf and we update global L and R structures. 

Observation 3 Every LHT generated by algorithm ff-fixup is an LST centered in a unique level 
graph. 

Proof. As described in (a) and (b) above, every LHT is generated using two triples which are 
shortest in the same level graph Tj. Moreover, since at least one of them must be centered in level 
i, the resulting LHT is an LST centered in level i. □ 

The last novelty in the algorithm is updating center arrays (Alg. [ 11 ] called at step [TO] - Alg. 
[IQI) in a similar way of ff-cleanup: Algorithm II II identifies the oldest level j related to the triples 
contained in Sb (Step [ 2 ]- Alg. [TT]l . If j > k then we reshape the distribution for 7 ' similarly to ff- 
CLEANUP (Steps [4] to [ 6 ]- Alg. Illh . Otherwise 7 ' is completely contained in level k and no reshaping 
is required (Step [ 8 ]- Alg. [TT]l . 


6 New Features in Algorithm ffd 


In this section, we discuss two challenges that arise when we attempt to generalize the level graph 
method used in Thorup (for APSP with unique SPs) to a fully dynamic APASP algorithm. Both 
are addressed by the algorithms in Section [5] as noted below. 

The bit /3 feature. The control-bit j3 was introduced (and only briefly described) in PR to avoid 
the processing of untouched historical triples. Here, we elaborate on this technique in more detail 
than PR, and we also describe how it helps in the more complex setting of the level tuple-system. 


Consider Figure [T] The ST 7 = {{xa, by),wt, count) is created in level k (Fig. 1(a)). At time{k), 
we have 7 G P* and also 7 € P with /3{'y) = 1. In a more recent level j < k, a shorter triple 
7 ' = {{xv, vy),wt', count'), with weight wt' < wt, that goes trough an updated vertex v is generated 


(Fig. 1(b)). Thus at time{j), we have 7 ' G P* and also 7 ' G P with / 3 ( 7 ') = 1; but 7 still appears in 


both P* and P as a historical triple. Finally, a new LST 7 " = {{xa', b'y),wt, count"), with the same 
weight as 7 , is generated in level i < j (Fig. 1(c)). Note that 7 " is only in P with I3{'y") = 0 and 
not in P*, as is the case of every LST that is not an ST. When an increase-only update removes v 
and the triple j', the algorithm needs to restore all the triples with shortest weight wt. But while 
7 is historical and does not require any additional extension, 7 " is only present in P and needs 
to be processed. Our ffd algorithm achieve this by checking the bit /3 associated to each of these 
triples. The algorithm will extract and process all the triples with /3 = 0 from P{x,y) (see StepISj 
Alg. [9] where triples are extracted using /3 = 0 as part of the key). These guarantees that a triple 
only present in P, or present in P and P* with different counts is never missed by the algorithm. 
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Fig. 1. The bit (3 feature 


The partial extension problem (PEP). Consider the update sequence described below and 
illnstrated Fig. [2j Here the STs 7 = {{xa, by), wt, count) and 7 = {{xa,cy),wt, count') are created 
in level k (Fig. 2(a)). Later, a left-extension to x' generates the STs 7 ' = {{x'x, by),wt', count) and 
y = {{x'x,cy),wt', count') in level j < k (Fig. 2(b)). Note that 7 , 7 , 7 ' and y are all present in 
P* and P at tinie{j). In a more recent level i < j, a decrease-only update on v generates a shorter 
triple 7 s = {{xv, vy),wts, counts) from x to y, with wts < wt going through v. In the same level, the 
triple 7 s is also extended to x' generating a triple 7 " = {{x'x,vy),wt", counts) shorter than 7 ' and 
7 ' (Fig. 2(c)). Thns at time(i), 7 , 7 , 7 ' and y remain in P* as historical triples. Then, in level h < i, 
an update on x” inserts the edges {x", x) and {x", c). This update generates an ST 7 '" = {{x"c, cy)) 
(shorter than {x"x,vy)) and also inserts x" G LC^{x,b)) since {x",x) is on a shortest path from x" 
to b; bnt it shonld not generate the triple of the form {x"x, by) because b is not on a shortest path 
from X to y at tinie{h) (Fig. 2(d)). 


X 

p 



y 


(a) level k 



y 

(b) level j < k 



y 


(c) level i < j 



Fig. 2. PEP instance (only centered STs are kept in each level) - all edge weights are unitary 


When an increase-only npdate removes v and the shortest triple 7 ^ from x to y, the algorithm 
needs to restore all historical triples with shortest weights from x to y. When 7 and 7 are restored, 
we need to perform suitable left extensions as follows. An extension to x" is needed only for 7 : in 
fact 7 should not be extended to x" becanse the £-tuple of the newly generated tuple is not an ST 
in the graph. On the other hand, no extension to x' is needed since both 7 ' and y will be restored 
(from HT to ST). Our algorithm needs to distinguish all of these cases correctly and efficiently. 
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In order to maintain both correctness and efficiency in this scenario for APASP, we use two new 
data structures (described in Section [4.2|] : (1) the historical distance matrix DL that allow us to 
efficiently determine the most recent level graph in which an HT was an ST (see for example Steps [5] 
and EH Alg. llOp . and (2) the HE sets LN and RN that allow us to efficiently identify exactly those 
new extensions that need to be performed (see for example Steps [Tl]~ ESI Alg. [9]). The methodology 
of these data structures is fully discussed in the description of ff-fixup (Section I5.2jl . Note that, 
the PEP doesn’t arise in Thorup because of the unique SP assumption: in fact when only a single 
SP of a given length is present in the graph for each pair of nodes, the algorithm can check for 
all the O(n^) paths maintained in each level and decide which one should be extended. Given the 
presence of multiple SPs in our setting, we cannot afford to look at each tuple in the tuple-system 


7 Correctness and Complexity 

In this section we will hrst prove the complexity bounds of our ff-update algorithm, then we will 
establish correctness. 

7.1 Complexity 

The complexity analysis of algorithm ffd is similar to that for the PR algorithm. We highlight the 
following new elements: 

1. As noted in Section 15.21 (see Observation [3]) , every triple created by ff-fixup is an LST in the 
level graph (PDG) in which is centered, and by the decremental only properties of level graphs, 
it will continue to be an LST in that level graph until it is removed. In contrast, PR can create 
LHTs by combining HTs not centered in any PDG. This results in an additional ©(log re) factor 
in the amortized bound there. 

2. We can bound the number of LHTs that contain a given vertex u as 0{z' ■ where z' is the 
number of active level graphs that contain vertex u and tuples passing through u (by Corollary 
E]). Given our level tuple-system, z' is clearly ©(logre). In PR, this bound is (z + z'^) where z is 
the number of active PDGs, and z' is the number of PDGs that contain v. 

3. We can show that the number of accesses to RN and LN^ outside of the newly created tuples, 
is worst-case 0{n ■ v*) per call to ff-update. The overhead given by the level data structures 
is ©(logre) for each access (see Lemma [5]) ■ These structures are not used in PR. 

Lemma 1. Let G be a graph after a sequence of calls to ff-update. Let z be the number of active 
level graphs (PDGs), and let z' < z be the number of level graphs that contain a given vertex v. 
Suppose that every HT in the tuple-system is an ST in some level graphs, and every LHT is an 
LST in some level graph. If re and m bound the number of vertices and edges, respectively, in any 
of these graphs, and if v* bounds the maximum number of distinct edges that lie on shortest paths 
through any given vertex in any of the these graphs, then: 

1. The number of LHTs in G’s tuple-system is at most 0{z ■ m ■ v*). 

2. The number of LHTs that contain a vertex v in G is 0{z' ■ re*^). 

Proof. Eor part 1, we bound the number of LHTs {xa,by) (across all weights) that can exist in G. 
The edge (x, a) can be chosen in m ways, and once we fix (x, a), the r-tuple (a, by) must be an ST 
in one of the Pj. Since (6, y) must lie on a shortest path through a centered in a graph Pi, that 
contains the r-tuple {a, by) of shortest weight in T), the number of different choices for {b,y) that 
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will then uniquely determine the tuple {xa,by), together with its weight, is z-iy*. Hence the number 
of LHTs in G’s tuple-system is 0{z ■ m ■ v*). 

For part 2, the number of LHTs that contain v as an internal vertex is simply the number of 
LSTs across the z' graphs that contains u, and this is 0(^z' ■ u*‘^). We now bound the number of 
LHTs (uo, by). There are n — 1 choices for the edge {v, a) and z' ■ u* choices for the r-tuple (a, by), 
hence the total number of such tuples is 0{z' ■ n ■ v*). The same bound holds for LHTs of the form 
{xa,bv). Since v* = I7(n), the result in part 2 follows. □ 

Corollary 1. At a given time step, let B be the maximum number of tuples in the tuple-system 
containing a path through a given vertex in a given level graph. Then, B = 0{v*‘^). 

Lemma 2. (a) - The cost for an ff-cleanup call on a node v when z' active levels contain triples 
through v is 0{z' ■ u*'^ ■ logn). 

(b) - The cost for a real ff-cleanup call is 0(p*^ • log^ n) 

(c) - The cost for a dummy ff-cleanup call is ■ logn). 

Proof, (a) - Algorithm 0] extracts all the LHTs that go through the update vertex from the heap 
He. Since the number of these LHTs is bounded by B at each level (by Corollary [T]) , the total cost 
is 0{z' ■ B\ogn) where z' is the number of active levels that contain triples through v. Algorithm 
El requires only O(logn) time for each step, except in step |4] where the cost is O(logn) for each 
triple extracted from He that goes through the updated vertex. Since the number of such triples 
is bounded by 0{z' ■ (by Lemma [T]), the worst-case cost for a call to Algorithm [6] within an 

FF-CLEANUP phase is 0{z' ■ • logn). In Algorithm El a triple can be added to heap He, or 

searched and removed from a constant number of priority queues among z! different active levels. 
Moreover, for the structures DL, RN and LN each triple spends a constant time to be unlinked 
and eventually to update the structures. Since, priority queue operations have a O(logn) cost and 
the number of triples examined is bounded by 0{z' ■ n*^), the complexity for Algorithm El is at 
most 0{z' ■ • logn). Thus an ff-cleanup call that operates on z' active levels requires at most 

0{z' ■ ■ logn). 

(b) - Since 2 ; < log 2 n, the cost for a real ff-cleanup call is 0{v*^ ■ log^ n) (by part (a)). 

(c) - For a dummy cleanup on a vertex w, ff-cleanup only needs to clean the local data structures 
in level center{w), where w is centered, and in the current level graph. In fact, let t be the current 
update step; in the dummy cleanup phase, we start with the node u that was updated at time t — 1 
(the most recent update before the current one). The node u received an update in the previous 
phase, thus it disappeared from all the levels older than level{t — 1) and, with it, all the LSTs 
containing u in these levels. Hence, all the triples containing u in the tuple-system must be LSTs 
in level{t — 1). We have at most B of them and ff-cleanup spends 0{B ■ logn) (considering the 
access to the data structures) to remove them. Then, the dummy update reinserts u only in the 
current graph. The next phase moves on the node u' updated at time t — 2. Again, all the tuples 
containing u' must be LSTs in level{t — 2) and eventually the current graph if they were inserted 
because of the previous dummy update on u. 

Suppose in fact that there is a tuple 7 that contains u' in another level (except the current 
graph). The tuple 7 cannot be in a level older than level{t — 2) because when u' was updated at 
time {t — 2), the cleanup algorithm removed all the tuples containing u' from any level older than 
t — 2. Moreover, a tuple containing u' present in a level younger than level{t — 2) could appear 
if and only if it was generated by any update more recent of t — 2 (in this case only the dummy 
update on u performed in the current graph). Thus a contradiction. 

This argument can be recursively applied to every other node in the sequence: in fact for the 
node u” updated at time {t — i) all the nodes updated in the interval [t — i-\-l,t — l] will be already 
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processed by ff-cleanup, leaving all the tuples containing u" only in level{t — i) and t. It follows 
that, for a dummy update, z' = 2. Thus the cost for a dummy ff-cleanup call is 0 (p*^ • logn) 
(by part (a)). □ 

Lemma 3. The. cost for a dummy ff-fixup call on a node v is • logn). 

Proof. Consider a dummy FF-FixuP applied to a vertex v in J\f. We only need to bound the cost 
for accessing the entries in the P*{x,y) and the cost of re-adding LSTs containing v, previously 
removed by the dummy ff-cleanup but still in the current graph after the dummy update. In 
fact the vertex v is removed by an earlier dummy ff-cleanup, and while this removes all the 
HPs containing the vertex v, it does not change any LST centered in any Pj that does not contain 
V. Hence these other LSTs will be present in the tuple-system with unchanged weight and count, 
when dummy FF-FIXUP is applied to v. Since for any pair x,y, the SP distance will not change 
after the dummy update, the dummy ff-fixup will only insert in the set S triples containing the 
node V for additional extension (this is accomplished by the check at Step [2] - Alg. [H followed 
by Steps [T8II221 Alg. [9]). Hence, only the LSTs containing v in the current level{t) graph will be 
processed and added to the tuple-system, and there are at most B of them (by Corollary [T]) . Thus 
a dummy FF-FIXUP for any v needs to access P* for each pair of nodes, and reinsert at most B 
tuples (containing v) in the current graph. Hence the overall complexity for a dummy ff-fixup is 
0 ((re^-|-H) • logn) = 0 (u*^ • logn). □ 

We now address the complexity of a real ff-fixup call. We first define the concept of a triple 
pair that will be used in lemma [5] to establish the bound for a real ff-fixup call. Finally, we 
complete our analysis by presenting a proof of Theorem [TJ 

Definition 1. If C^[i\ > 1 then ( 7 ,i) is a triple pair in the tuple-system. If ( 7 ,i) is not a triple 
pair in the tuple-system at the start of step t but is a triple pair after the update at time step t, 
then ( 7 , i) is a newly created triple pair at time step t. 

Lemma 4. At a given time step, let D be the number of triple pairs in the level tuple-system. Then, 

1. The value of D is at most 0{m ■ v* ■ logn). 

2. The space used is 0{m ■ u* ■ logn). 

Proof. 1. Every > 1 represents a distinct LST in Pi, hence the result follows since the number 
of levels is O(logn) and the number of LSTs in a graph is 0(n* • m*). 

2. Since every triple is of size 0(1), the memory used by our ffd algorithm is dominated by D, 
and result follows from 1 . □ 

Lemma 5. The cost for a real ff-fixup call is 0{u*^ ■ log^ n -|- X • log n) , where X is the number 
of newly created triple pairs after the update step. 

Proof. A triple is accessed only a constant number of time during ff-fixup for a cost of O(logn), 
so it suffices to establish that the number of existing triples accessed during the call is 0 (n*^ -log n). 

There are only O(n^) accesses to triples in the call to Algorithm 8 in line 2 of ff-fixup since 
O(n^) entries in the global P*{x,y) structures are accessed to initialize Hf. This takes 0(n^ Tog re) 
time after considering the O(logre) cost per data structure operation. We now address the accesses 
made in the main loop. We will distinguish two cases and they will be charged to X as follows. 

1: fd{j) = 0 - In the main loop of Algorithm 7, starting in Step 3, any triple 7 that is accessed with 
/ 3 ( 7 ) = 0 is an LST at some level i where it is not identified as an ST in T). During this call, in 
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Steps [3HI21 Alg. [9] (or Steps [T8H221 Alg. [9] if the distance for the endpoints of 7 did not change) 7 
is added as an ST in level i, and will never be removed as an ST for level i nntil it is removed from 
the tuple system (due to the fact that Fi is a purely decremental graph). Since 7 with / 3 ( 7 ) = 0 is 
a newly added triple to level i, then the pair ( 7 , i) is a newly created triple pair at step t. Hence, 
we can charge ( 7 , i) to X in this call of ff-fixup. 

2 : / 3 ( 7 ) = 1 - We now consider triples accessed that have /? = 1. This is the most nontrivial part of 
our analysis since even though any such triple 7 must exist with the same count in every level in 
both P and P*, we may still need to form some extensions since the triple may have been an HT 
when extension vertices were updated, and hence these extension may not have been performed. 
Here is where the LN and RN sets are accessed, and we now analyze the cost of these accesses. 
(The correctness of the associated steps is analyzed in the next section.) 

Let j be the most recent level in which 7 was an ST in G and assume we are dealing with left 
extensions (right extensions are symmetrical). Now that 7 is restored, the only case (Steps [TT] to 
[THl Alg. [ 9 ]) in which we need to process it is when there exists a left extension for the Atuple of 7 
to a node x' centered in a level i more recent than j. In fact, the LST generated by this extension 
will appear for the hrst time centered in level i, hence the pair ( 7 , i) is a newly created triple pair 
at step t and we can charge its creation to X. We now show how our HE sets efficiently handle this 
case. By steps [H] to [161 Alg. [H our algorithm only processes a restored triple 7 with / 3 ( 7 ) = 1 when 
it has at least one centered extension in some active level younger than the level in which 7 was 
shortest for the last time. We can bound the total computation for these steps as follows: for a given 
X, RN{x, y, wt) contains a node b for every incoming edge to y in one of the SSSP dags (historical 
and shortest) rooted at x. Since we can extend in at most O(logn) active levels during any update 
and the size of a single dag is at most v*, these steps take time 0{v* -nlogn) throughout the entire 
update computation. □ 

We can now establish the proof of our main theorem. 

Proof of Theorem [T1 Consider a sequence X oi r = Q{n) calls to algorithm ffd. Recall that 
the data structure is reconstructed after every 2n steps, so we can assume r = 0{n). These r 
calls to FFD make r real calls to ff-update, and also make additional dummy updates. As in PR, 
across the r real updates in A, the algorithm performs 0(r logn) dummy updates. This is because 
r/ 2 ^ real updates are performed at level k during the entire computation, and each such update is 
accompanied by 2 ^ — 1 dummy updates. So, across all real updates there are 0{r) dummy updates 
per level, adding up to 0 (r log n) in total, across the O(logn) levels. 

When FF-CLEANUP is called on a vertex v for a dummy update, z' = 2 since v can be present 
only in the most recent current level and the level at which it is centered. (This is because every 
vertex that was centered at a more recent level than v has already been subjected to a dummy 
update, and hence all of these vertices are now centered in the current level.) Thus, by Lemma O 
each FF-CLEANUP for a dummy update has cost 0{B -logn). By Lemma[3l a call to ff-fixup for a 
dummy update has cost 0{v*^-log n). Thus the total cost is 0{{u*^ -log n) -r log n) across all dummy 
updates. Also, the number of tuples accessed by all of the dummy update calls to ff-cleanup, 
and hence the number of tuples removed by all dummy updates, is 0{r ■ • logn). 

For the real calls to ff-fixup, let Xi be the number of newly added triple pairs in the ith real call 
to FF-FIXUP. Then by LemmaO the cost of this iih. call is 0(n*^-log^ n-|-W-log n). Let X = Yll=i 
Hence the total cost for the r real calls to ff-fixup is 0(r • Tog^ n-|-A-log n). We now bound X 
as follows: X is no more than the maximum number of triples that can remain in the system after X 
is executed, plus the number of tuples Y removed from the tuple-system. Tuples are removed only 
in calls to ff-cleanup. The total number removed by rlogn dummy calls is 0(r • logn • n*^) (by 
Lemma[2]). The total number removed by the r real calls is 0{r ■ • logn) (by Lemma[2]). Hence 
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V = 0(r ■ u*'^ ■ logn). Clearly the maximum number of triples in the tuple-system is no more than 
D, which counts the number of triple pairs, and we have D = 0{m ■ v* ■ logn) = 0{v? ■ v* ■ logn) 
(by LemmaH]). Since r = 0(n), we have D = 0(r ■ n ■ u* ■ logn), and this is dominated by Y since 
V* = Q{n). Hence the cost of the r calls to ffd is 0{r ■ • log^ n) (after factoring in the O(logn) 

cost per tuple access), and hence the amortized cost of each call to ff-update is 0(n*^ • log^ n). 

7.2 Correctness 

For the correctness, we assume that all the global and local data structures are correct before the 
update, and we will show the correctness of them after the update. 

Correctness of Cleannp - The correctness of ff-cleanup is established in Lemma [H We will 
prove that all paths containing the updated vertex v are removed from the tuple-system. Moreover, 
the center of each triple is restored, if necessary, to the level containing the most recently updated 
node on any path in this triple. Note that (as in |8ll9l21j i at the end of the cleanup phase, the 
global structures P and P* may not have all the LHTs in G \ {u}. 

Lemma 6. At the end of the cleanup phase triggered by an update on a vertex v, every LHP that 
goes through v is removed from the global structures. Moreover, in each level graph Pi, each SP that 
goes through v is removed from P*. For each level i, the local structures L*, R*, RC* and LC* 
contain the eorrect extensions; the global struetures L and R eontain the eorrect extensions, for 
each r-tuple and £-tuple respeetively, and the struetures RN and LN eontain only nodes assoeiated 
with tuples in P. The DL structure only contains historical distances represented by at least one 
path in the updated graph. Finally, every triple in P and P* has the correct updated center for the 
graph G \ {u}. 

Proof. To prove the lemma statement, we use a loop invariant on the while loop in Step [3] of 
Algorithmic We show that the while loop maintains the following invariants. 

Loop Invariant: At the start of each iteration of the while loop in Step [3] of Algorithmic assume 
that the hrst triple to be extracted from He and processed has min-key = \wt,x,y]. Then the 
following properties hold about the tuple-system and He. 

1. For any a, 6 G H, if G contains Cab LHPs of weight wt of the form {xa,by) passing through v, 
then He contains a triple 7 = {{xa,by),wt,Cah) with key [wt,x,y\ already processed: the Cab 
LHPs through v are not present in the tuple-system. 

2. Let [wt,x,y\ be the last key extracted from He and processed before \wt,x,y]. For any key 

\wti,xi,yi] < \wt,x,y\, let G contain c > 0 number of LHPs of weight wti of the form 
(xix, 6 i?/i). Further, let c„ (resp. Ca) denote the number of such LHPs that pass through 
V (resp. do not pass through v). Here = c. For every extension x' G L{xi,biyi), let 

wt' = wti -|-w(x',xi) be the weight of the extended triple fx'xi,biyi). Then, (the following 
assertions are similar for y' G R{xiai,yi)) 

Global Data Structures: 

(a) if c > there is a triple in P{x',yi) of the form {x'xi,biyi) and weight wt' representing 
c — Cy LHPs. Moreover, its center is updated according to the last update on any path 
represented by the triple. If c = c.;, there is no such triple in P{x',yi). 

(b) If a triple of the form {x'xi,biyi) and weight wt' is present as an HT in P*{x',yi), then it 
represents the exact same number of LHPs c — Cy of the corresponding triple in P(x',yi). 
This is exactly the number of HPs of the form (x'xi, biyi) and weight wt' in G \ {u}. 

(c) x' G L{xi,biyi), yi G R{x'xi,bi), and {x'xi,biyi) G Marked-Tuples iff Cy > 0. 
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(d) A triple corresponding to {x'xi,hiyi) with weight wt' and counts is in He- A similar 
assertion holds for y' G i?(xiai,yi). 

(e) The structure RN{x',yi,wt') contains a node b iff at least one path of the form {x'x,byi) 
and weight wt' is still represented by a triple in P{x',yi). A similar assertion holds for a 
node a in LN{x',yi,wt'). 

(f) If there is no HT of the form {x'x, biyi) and weight wt' in P*{x' , yi) then the entry DL{x' , yi) 
with weight wt' does not exists. 

Local Data Structures: for each level j, let Cj be the number of LSPs of the form {x'xi, biyi) 
and weight wt' centered in Pj and let Cj{v) be the ones that go through v. Thus c = J2j Cj and 
Cv = Yjj Cj(v)- Then, 

(g) the value of C.y[j], where 7 is the triple of the form {x'xi, biyi) and weight wt' in P{x', yi), 
is Cj — Cj{v). 

(h) If a triple 7 of the form {x'xi, biyi) and weight wt' is present as an HT in P*, then P*{x', yi) 
represents only Cj — Cj{v) paths. If Cj — Cj{v) = 0 then the link to 7 is removed from dietj. 
Moreover, x' G L*{xi^yi) (respectively LC*{xi^yi) if x' is centered in Pj) iff x' is part 
of a shortest path of the form {x'xi, xyi) centered in Pj. A similar statement holds for 
yi G R*{x',bi) (respectively RC*{x',bi) if yi is centered in Pj). 

3. For any key [wt 2 , X 2 ,y 2 ] > [wt,x,y], let G contain c > 0 LHPs of weight wt 2 of the form 
(x 2 a 2 , 622 / 2 )- Further, let c„ (resp. Cy) denote the number of such LHPs that pass through v 
(resp. do not pass through u). Here Cy + Cy = c. Then the tuple {x2a2, 622 / 2 ) G Marked-Tuples, 
iff Cci > 0 and a triple for {x 2 a 2 , 622 / 2 ) is present in Hy 

Initialization: We start by showing that the invariants hold before the first loop iteration. The 
min-key triple in Hy has key [0,u,u]. Invariant assertion [1] holds since we inserted into Hy the 
trivial triple of weight 0 corresponding to the vertex v and that is the only triple of such key. 
Moreover, since we do not represent trivial paths containing the single vertex, no counts need to be 
decremented. Since we assume positive edge weights, there are no LHPs in G of weight less than 
zero. Thus all the points of invariant assertion [2] hold trivially. Invariant assertion [3] holds since Hy 
does not contain any triple of weight > 0 and we initialized Marked-Tuples to empty. 
Maintenance: Assume that the invariants are true before an iteration k of the loop. We prove 
that the invariant assertions remain true before the next iteration k + 1. Let the min-key triple 
at the beginning of the fc-th iteration be [rct/j, 2 /fc]. By invariant assertion [H we know that for 
any a*, bj, if there exists a triple 7 of the form (x^Oj, bjyk) of weight wtk representing count paths 
containing v, then it is present in Hy. Now consider the set of triples with key [wtk,Xk, yk\ which we 
extract in the set S (Step 01 Algorithm 0]) . We consider left-extensions of triples in S; symmetric 
arguments apply for right-extensions. Consider for a particular 6 the set 5^ C S' of triples of the 
form (xfc—, 62 /A:), and let /count' denote the sum of the counts of the paths represented by triples in 
Sb- Let x' G L{xk, byk) be a left extension; our goal is to generate the triple 7' of the form {x'xk, byk) 
with count /count' and weight wt' = wtk + ’w(x',xa;), and an associated vector C{'y') that specifies 
the distribution of paths represented by j' level by level. These paths will be then removed by the 
algorithm. However, we generate such triple only if it has not been generated by a right-extension 
of another set of paths by checking the Marked-Tuples structure: we observe that the paths of the 
form {x'xk, byk) can be generated by right extending to yk the set of triples of the form {x'xk, x 6 ). 
Without loss of generality assume that the triples of the form {x'xk, x 6 ) have a key which is greater 
than the key [wtk, Xk,yk]- Thus, at the beginning of the fc-th iteration, by invariant assertion[3l we 
know that {x'xk,byk) ^ Marked-Tuples. StepO Alg. 0] creates a triple 7' of the form {x'xk,byk) of 
weight wt' and /count'. 

The set of triples in Sb can have different centers and we are going to remove (level by level) 
paths represented by yh To perform this task we consider the vector C.yir. it contains the full 
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distribution of the triple 7" G P{x'yk) of the form {{x'xk,byk),wt') and indicates the oldest level 
j in which 7 '' was generated for the first time. This level is exactly argmaxj(C'.y//[j] 7 ^ 0) and it 
is identified in Step [3] - Alg. [5j All the paths represented in Sb that are centered in some level m 
older or equal to j were extended for the first time in level j to generate y. Moreover, each path 
centered in a level i younger than j was extended in level i itself. Thus, we can compute a new 
center vector C^' (according to the distribution in Cy/) of the paths containing v that we want to 
delete at each active level, as in steps [5] to [ 8 ] - Alg. [ 6 l In step [9]- Alg. [ 6 ] the vector Cy is updated: 
the paths are removed level by level according to the new distribution. This establishes invariant 
assertion 

The triple 7' is immediately added to He with C^i for further extensions (Step. [7]- Alg. [5]). This 
establishes invariant assertions IMl Thus we reduce the counts of 7' in P{x',yk) by /count (Step. 
El- AlglS]) and we set the new center for the remaining tuple j" in P(x',yk) establishing invariant 
assertion [ 2 al Steps [10] to [13] - Alg. [5] check if there is any path of the form {x—,by) that can use 
x' as an extension. In this case we add 'j' to the Marked-Tuples. If not, we safely remove the left 
and right extension [x' and y) from the tuple-system. This establishes invariant assertion [23 If 7 ^ 
is an HT in P*{x',yk), we decrement its count (Step. [19]- Alg[5]) establishing invariant assertion 
[26] In steps fTTl to fT71 - Alg. [5] we use the double links between 6 G RN{x',yk,wt') and tuples 
to efficiently check if there are other triples linked to 6 ; if not we remove 6 from RN(x',yk,wt') 
establishing invariant assertion!^ Using a similar double link method with the structure DL{x', yk), 
we establish invariant assertion 2/ after step [21] - Alg. [5] 

To operate in the local data structures we require 7' to be an HT in P*{x',yk)- Using the 
previously created vector we reduce the count associated with 7' G lP*{x',yk) for each level i 
(Step [23] - Alg. [5]). After the above step, if there are no paths left in P*{x',yi) then there are no 
STs of the form {x'xi,byi) centered in level i. In this case we remove the extension x' and y^ from 
the local structures of level i. This is done in steps [24] to [34] -Alg. [5] in case x' is not centered in 
level i, then any path in 7' centered in level i is generated by a node centered in level i located 
between Xk and yk- Thus if any SP from x' to yk (that uses {x',Xk) as a first edge) remains in in 
Pi, it must be also counted in P*{xk,yk)- Thus, we remove x' from L*{xk,yk) only if P*{xk,yk) is 
empty. In the case x' is centered in level i and P*{xk,yk) is empty, x' could still be the extension of 
other paths from Xk to yk centered in levels older than i. The algorithm checks them all and if they 
do not exist in older levels we can safely remove x' from LC*{xk,yk) (Step [32l Alg. [5]). A similar 
argument holds for the right extension yi- This establishes invariant assertion 126.1 and completes 
claim [2 

When any triple is generated by a left extension (or symmetrically right extension), it is inserted 
into He as well as into Marked-Tuples. This establishes invariant assertion [3] at the beginning of 
the [k + l)-th iteration. 

Finally, to see that invariant assertion [T] holds at the beginning of the (/c-|-l)-th iteration, let the 
min-key at the {k + l)-th iteration be \wtkJei,Xk+i,yk+i[- Observe that triples with weight wtk+i 
starting with Xk+i and ending in yk+i can be created either by left extending or right extending 
the triples of smaller weight. And since for each of iteration < k, invariant assertion [2] holds for any 
extension, we conclude that invariant assertion [ 1 ] holds at the beginning of the {k + l)-th iteration. 
This concludes our maintenance step. 

Termination: The condition to exit the loop is He = 0. Because invariant assertion [1] maintains 
in He all the triples already processed, then He = 0 implies that there are no other triples to 
extend in the graph G that contain the updated node v. Moreover, because of invariant assertion 
[T] every triple containing the node v inserted into He has been correctly decremented from the 
tuple-system. Remaining triples have the correct update center because of invariant [2a] Finally, for 
invariant assertions and [Ml the structures L *, LC*, R*, RC* are correctly maintained for every 
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active level i and the paths are surgically removed only from the levels in which they are centered. 
This completes the proof. □ 

Correctness of Fixup - For the fixup phase, we need to show that the triples generated by our 
algorithm are sufficient to maintain all the ST and LST in the current graph G. As in PR, we first 
show in the following lemma that ff-fixup computes all the correct distances for each pair of nodes 
in the updated graph. Finally, we show that data structures and counts are correctly maintained 
at the end of the algorithm (Lemma [8]). 

Lemma 7. For every pair of nodes {x,y), let 7 = {{xa, by), wt, count) be one of the min-weight 
triples from x to y extracted from Hf during ff-fixup. Then wt is the shortest path distance from 
X to y in G after the update. 

Proof. Suppose that the lemma is violated. Thus, there will be an extraction from Hf during 
FF-FIXUP such that the set of extracted triples S', of weight wt is not shortest in G after the 
update. Consider the earliest of these events when S' is extracted from Hf. Since S' is not a set 
of STs from x to y, there is at least one shorter tuple from x and y in the updated graph. Let 
7' = {{xa', b'y),wt, count) be this triple that represents at least one shortest path from x to y, with 
wt < wt. Since S' is extracted from Hf before any other triple from x to y, 7 ' cannot he in Hf at 
any time during FF-FIXUP. Hence, it is also not present in P{x, y) as an LST at the beginning of 
the algorithm, otherwise it (or another triple with the same weight) would be placed in Hf hy step 
[ 2 ]- Alg. [3 Moreover, if 7' is a single edge (trivial triple), then it was already an LST in G present 
in P{x, y) before the update, and it is added to Hf hy step [10]- Alg. jS] moreover since all the edges 
incident to v are added to Hf during steps [2| to [3 of Alg. [ 8 l then 7' must represent SPs of at least 
two edges. We define left{y') as the set of LSTs of the form {{xa', Cib'),wt — w(5', y), counted) that 
represent all the LSPs in the left tuple {{xa',b'),wt — w{b',y))-, similarly we define right{y') as the 
set of LSTs of the form {{a'dj,b'y),wt — w{x, a'), county.) that represent all the LSPs in the right 
tuple {{a',b'y),wt — 'w{x,a')). 

Observe that since 7 ' is an ST, all the LSTs in left{y') and right{y') are also STs. A triple in 
left{'y') and a triple in right{'y') cannot be present in P* together at the beginning of ff-fixup. In 
fact, if at least one triple from both sets is present in P* at the beginning of ff-fixup, then the last 
one inserted during the fixup phase triggered during the previous update, would have generated an 
LST of the form {{xa',b'y),wt) automatically inserted, and thus present, in P at the beginning of 
the current fixup phase (a contradiction). Thus either there is no triple represented by left{y') in 
P*, or there is no triple represented by right{y') in P*. 

Assume w.l.o.g. that the set of triples in right{y') is placed into P* after left{'y') by ff-fixup. 
Since edge weights are positive, wt — w(x, a') < wt < wt, and because all the extractions before 7 
were correct, then the triples in right{y') were correctly extracted from Hf and placed in P* before 
the wrong extraction of S'. Let i be the level in which left{'y') is centered, and let j be the level in 
which right{y') is centered. By the assumptions, all the triples in left{y') are in P* and we need 
to distinguish 3 cases: 

1. if j = i, then ff-fixup generates the tuple {{xa',b'y),wt) in the same level and place it in P 
and Hf. 

2. if i > j, the algorithms ff-fixup extends the set right{'y') to all nodes in L*{a',b') for every 
i P j (see Steps 13 to dH - Alg. [T 0 ]l . Thus, since left{y') is centered in some level i > j, the node 
x is a valid extension in L*{a', b'), making the generated y' an LST in Pj that will be placed in 
P{x,y) and also into Hf (during Step fTOl - Alg. fTOl). 
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3. if j > i, then x was inserted in a level younger than i. In fact, all the paths from a' to b' must be 
the same in right{'y') and otherwise the center of right{'y') should be i. Hence, the only 

case when j > i is when the last update on left{'y') is on the node x in a level i younger than 
j. Thus X G LC*{a',b'). But ff-fixup extends right{'y') to all nodes in LC*{a',b') for every 
i < j, placing the generated LST 7 ' in P{x,y) and also into Hj (see Steps [IT] to [19] - Alg. [Toll . 

Thus the algorithm would generate the tuple {{xa', b'y), wt) (as a left extension) and place it 
in P and Hf (because all the triples in left^j') are already in P*). Therefore, in all cases, a tuple 
{{xa', b'y), wt) should have been extracted from before any triple in S'. A contradiction. □ 

Lemma 8 . After the execution o/ff-fixup, for any {x, y) G V, the sets P*{x, y) (P{x, y)) contains 
all the SPs (LSPs) from x to y in the updated graph. Also, the global structures L,R and the local 
structures P*, L*, R*, LC *, RC* and dicti for each level i are correctly maintained. The structures 
RN and LN are updated according to the newly identified tuples. The DL structure contains the 
updated distance for each pair of nodes in the current graph. Finally, the center of each new triple 
is updated. 

Proof. We prove the lemma statement by showing the following loop invariant. Let G' be the graph 
after the update. 

Loop Invariant: At the start of each iteration of the while loop in Step [3] of ff-fixup, assume 
that the first triple in T/j to be extracted and processed has min-key = [wt, x, y]. Then the following 
properties hold about the tuple-system and Hf. 

1. For any a,b £ V, if G' contains Cab SPs of form {xa,by) and weight wt, then Hf contains 
a triple of form {xa, by) and weight wt to be extracted and processed. Further, a triple 7 = 
{{xa,by),wt,Cab) is present in P{x,y). 

2. Let [wt,x,y] be the last key extracted from Hf and processed before [wt,x,y]. For any key 
[wti, xi,yi] < [wt,x,y], let G' contain c > 0 number of LHPs of weight wti of the form 
{xiai,biyi). Further, let Cnew (resp. Coid) denote the number of these LHPs that are new (resp. 
not new). Here Cmw + Coid = c. If Cnew > 0 then, 

Global Data Structures: 

(a) there is an LHT 7 in P{xi,yi) of the form {xiai,biyi) and weight wti that represents c 
LHPs, with an updated center defined by the last update on any of the paths represented 
by the LHT. 

(b) If a triple of the form (xioi, biyi) and weight wti is present as an HT in P*, then it represents 
the exact same count of c HPs of its corresponding triple in P. This is exactly the number 
of HPs of the form {xiai, biyi) and weight wti in G'. Its control bit j3 is set to 1. 

(c) xi G L{ai,hiyi), yi G R{xiai,bi). Further, {xiai,biyi) G Marked-Tuples iff Coid > 0. 

(d) If (d{'f) = 0 or j3{'f) = 1 and there is an extension x' G L*{xi,yi) that generates a centered 
LST in a level j, an LHT corresponding to {x'xi, biyi) with weight wt' = wti + w{x', xi) > 
wt and counts equal to the sum of new paths represented by its constituents, is in Hf and 
P. A similar assertion holds for an extension y' G R*{xi,yi). 

(e) The structure RN{xi,yi,wti) contains a node b iff at least one path of the form {xi x, byi) 
and weight wti is represented by a triple in P{xi,yi). A similar assertion holds for a node 
o in LN{xi,yi,wti). 

(f) The entry DL{xi,yi) with weight wti is updated to the current level. 

Local Data Structures: for each level j, let Cj be the number of SPs of the form (xioi, biyi) 
and weight wti centered in Pj and let Cj{n) be the new ones discovered bythe algorithm. Thus 
c = Ei Cj and C-new — Yfj Cj{n). Then, 
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(g) the value of C^\j], where 7 is the triple of the form (xiui, 6iyi) and weight wti in P(xi, 7/1), 
is Cj. 

(h) If a triple 7 of the form (xioi, hiyi) and weight wti is present as an HT in P*, then P*{xl, yl) 
represents Cj paths. A link to 7 in P is present in dictj. Moreover, xi € L*{ai,yi) (respec¬ 
tively LC*{ai,yi) if xi is centered in Pj). A similar statement holds for yi G R*{xi^bi) 
(respectively RC*{xi,hi) if yi is centered in Pj). 

3 . For any key [wt2,X2,y2] > [wt,x,y], let G' contain c > 0 number of LHPs of weight wt2 of the 
form (x2a2, 627/2)- Further, let Cmw (resp. Coid) denote the number of such LHPs that are new 
(resp. not new). Here Cnew + Coid = c. Then the tuple {x2a2, 622/2) £ Marked-Tuples, iff Coid > 0 
and Cnew paths have been added to P/ by some earlier iteration of the while loop. 

Initialization and Maintenance for the invariant assertions above are similar to the proof of 
Lemma [H 

Termination: The condition to exit the loop is = 0. Because invariant assertion [1] maintains 
in iLj the first triple to be extracted and processed, then iLj = 0 implies that there are no triples, 
formed by a valid left or right extension, that contain new SPs or LSPs, that need to be added or 
restored in the graph G. Moreover, because of invariant assertions [ 2 ^ and [13 every triple containing 
the node v, extracted and processed before Pfj = has been added or restored with its correct count 
in the tuple-system. Finally, for invariant assertions ( 2 ^ and \ 2 hl the sets L,R and L*, LG*, R*, RG* 
for each level, are correctly maintained. This completes the proof of the loop invariant. 

By Lemma m all the SP distances in G' are placed in and processed by the algorithm. 
Hence, after Algorithm [ 7 | is executed, every SP in G' is in its corresponding P* by the invariant 
of Lemma [8l Since every LST of the form (xa, by) in G' is formed by a left extension of a set of 
STs of the form {ax, by) (Step[ 7 |- Algorithm [ 7 ]) , or a right extension of a set of the form {xa, xb) 
(analogous steps for right extensions), and all the STs are correctly maintained and extendend (by 
the invariant of LemmajS]), then all the LSTs are correctly maintained at the end of ff-fixup. This 
completes the proof of the Lemma. □ 
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